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We present an exact condition for the existence of gapless quasiparticle excitations in NJL models 
of color superconducting quark matter with a quark-quark interaction in the scalar color-antitriplet 
channel. The condition can be represented by a rotated ellipse in the plane of mass and chemical 
potential differences for the paired quark fields. 
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At high baryon density and low temperature mat¬ 
ter is believed to be in a color superconducting state, 
which is characterized by condensates of quark Cooper 
pairs Hiii A superconducting phase typically has 
an energy gap in the density of states, which corresponds 
to the lowest excitation energy of a quasiparticle pair. 
However, if the difference between the Fermi momenta of 
the paired quarks is sufficiently lar^ gapless quasipar¬ 
ticle excitations could exist S S 0 ■ The presence of 
gapless phases could have observable consequences, e.g., 
the high specific heat and neutrino emissivitv could af¬ 
fect the cooling behavior of compact stars [2|. It has 
been found, however, that gapless phases might suffer 
from a chromomagnetic instability [If] . UlL ll^ 

and it is presently unclear whether gapless phases ap¬ 
pear at temperatures relevant for compact star evolu¬ 
tion 0 [l^ [l^ I 12 I . It is therefore important to improve 
the understanding of gapless phases. In this paper we de¬ 
rive an exact condition for the existence of gapless excita¬ 
tions in the frequently used Nambu-Jona-Lasinio (NJL) 
model of color superconducting quark matter. A quali¬ 
tatively useful graphical representation of the condition, 
and some well-known approximations are also presented. 

The most dense environment where quark matter is 
expected to exist is in the core of neutron stars, which 
are subject to a gravitational instability that limits the 
maximum density to ~ 10^® g/cm^ [ 13 . This corre¬ 
sponds to a maximum quark number chemical potential 
of /i ^ 500 MeV and a maximum baryon number density 
oiriB ~ 10 no, where no = 0.17 fm“^ is the baryon num¬ 
ber density in nuclear matter. Since the charm quark 
mass is higher than the maximum chemical potential, it 
is sufficient to consider up (it), down (d) and strange (s) 
quarks. The quark spinors are 

Q — {gpur-, '4^ugi ‘4^ub: '4^dr^ '^dgi '4^dbt '4^ sr •, 4^sgt 4^sh') : ( 1 ) 

where r, g and b represent red, green and blue colors. The 
NJL model of superconducting quark matter is based on 
effective point-like four-fermion interactions and is de¬ 
scribed in, e.g., [HElIIlEi Here we repeat some of 


the essential points. The Lagrangian density is 

Ceff = q{i^ - m + fL'J°)q + Cqq + Cqq, (2) 

where m = diag/(TO„, mq, rris) is the current quark mass 
matrix in flavor space. Lqq and Lqq are the effective 
interaction terms, which are used at mean-field level in 
Hartree approximation. Explicitely 
8 

C-qq = Gs'^ \{qTaqf + {qilbTaqf 

a—0 

- -fs^[det/(g(l-k 75)9)-kdet/(g(l-75)9)], ( 3 ) 
Cqq = GD^{qh5TahCq'^){q'^Ci"/5TaXbq), (4) 

a,b= 2 , 5,7 


where Tq and A;, are the antisymmetric Gell-Mann matri¬ 
ces acting in, respectively, flavor and color space. Gs, K 
and Gd are coupling constants that must be determined 
by experiments. 

The quark-quark interaction term, Cqq , gives rise to su¬ 
perconducting condensates, Sab = (q'^CjbTaXbq), which 
break SU(3)c and U(l) symmetry. The symmetries of 
C corresponds to a conserved chromo-electromagnetic 
charge. The associated chemical potential is 
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+ T 3 X 3 + TsXs. (5) 


Here, fi is the quark-number chemical potential, /ig the 
positive electric-charge chemical potential, while /13 and 
/ig are color-charge chemical potentials. By linearizing 
@ in the quark-quark (diquark) gaps, Aab = ‘^GoSab, 
and the quark-antiquark (chiral) gaps, cfi = —AGs {qiqt), 
a grand canonical thermodynamic potential can be ob¬ 


tained by standard methods |l,5llHlll7l ll 
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Here, En{p, p; fiQ, ^. 3 ,^ 8 , 4'u,(t>d,(t>s, ^ud •) ^US •) Ads) are 
the quasiparticle dispersion relations, ^lep is the contri¬ 
bution from leptons (e.g., electrons, muons and the cor¬ 
responding neutrino flavors) and Hq is the vacuum, i.e., 
H(0, 0) = 0. It should be noted that ® is an even func¬ 
tion of En, so the signs of the dispersion relations are 
arbitrary. We therefore follow the standard convention 
that all states below the Fermi surface {En < 0) are oc¬ 
cupied, and only positive energy states are considered. 
In jnj the diquark gaps are denoted with flavor indices. 
This can readily be done by considering the color and 
flavor structure of the Gell-Mann matrices 


Fermi momenta of these three species are approximately 
equal , and no such gap less modes have been found in 
numerical evaluations jl^ [la, lnl . A proof has turned 
out to be difflcult to obtain due to the complexity of the 
characteristic polynomial of the 12x12 matrix. We there¬ 
fore leave this analysis to a future publication. Here the 
4x4 matrices are considered. The characteristic polyno¬ 
mials of these matrices can be written as 

En + O'SEn + 0‘2En + 0,iEn + Oq . (11) 

The ag coefficient of the polynomial is 113 
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and Aab = 0 if a yf 6 0 . The chiral gaps and the 
diquark gaps are variational parameters that are deter¬ 
mined by minimization of ®. The constituent quark 
masses are 

K 

Mi = mt + (l)i + ( 10 ) 

where (i, j, k) is any permutation of (u, d, s). 

In QCD, a color superconducting ground state is au¬ 
tomatically color neutral due to the generation of gluon 
condensates in one or more of the eight components of 
the gluon field. In NJL models there are no gauge fields 
that neutralize the color charge dynamically, because the 
gluons have been replaced by effective point-like quark- 
antiquark m and quark-quark 0 interactions. Color 
neutrality must therefore be enforced by solving for the 
charge chemical potentials, fiq, ps and /ig, such that 
the corresponding charge densities, Ua = (ip^Taip) = 
—dVl/dpa, are zero | 2 n| . 

The values of the gaps and the (charge) chemical po¬ 
tentials depend on the coupling constants (Gg, K and 
Gd), the current quark masses {niu, md and rus) and 
the regularization method. These input parameters are 
fitted to low-density hadronic results and are therefore 
only approximately known. In addition, approximations 
are frequently used in order to simplify the evaluation 
of ®. In this context it would be useful to have a math¬ 
ematically exact condition for the apperance of gapless 
quasiparticle dispersion relations, without reference to 
specific input parameters and further assumptions. This 
condition is presented below. 

The dispersion relations, if„, are eigenvalues of six 4x4 
matrices and one 12x12 matrix [lallflllS' Disregarding 
the signs, three 4x4 matrices and six of the twelve eigen¬ 
values of the 12x12 matrix remain (3 x 4 -|- 6 = 18). The 
12 x 12 matrix corresponds to ur — dg — sb pairing and 
the three 4x4 matrices correspond to ug — dr, ub — sr 
and db — sg pairing. There are strong indications that 
the ur — dg — sb modes are never gapless, because the 


ao=p‘^+ [mI + M] + 2A}j - 
+ {piapjp + MiMj + Afj + PiaMj + PjgMi) 

X {piaPjP + M^Mj + Afj - PiaMj - PjfdMi) , (12) 

for quark flavors {i, j) and colors (a, /3). The chemical 
potential, pia, for a quark held with flavor i and color a 
can be extracted from ® {p is diagonal in color and fla¬ 
vor space). A gapless dispersion relation is characterized 
by En{p) = 0 for some real value(s) of p when A^- ^ 0. 
This requires that ag{p) = 0 has at least one real root. 
The solutions are 


p^ = p^ + 6p^ -M^- SM^ - A^ 
±2^{pSp - MSM)^ - A 2(^2 _ sM2y 

Here we have introduced the quantities 


M = {M, + Mj)/2, SM = {Mi - Mj)/2, (14) 
p = {pia T — {Pia (^^) 

and A = A^. The indices in m can be omitted without 
ambiguity, since we are dealing with two-species pairing. 
Observe that the masses and chemical potentials of the 
paired quark fields are M±SM and p±Sp. A real square 
root in la requires that 


A < A9 


\p6p — M5M\ 
^Jp? - SM^ ’ 


(16) 


and a positive solution for p^ requires that 


-t SM^ + A^-p^- Sp^ 

< 2sJ{pSp - M6M)^ - A^{p^ - SM^). (17) 


Inequality can be represented with a rotated el¬ 
lipse in the SM — 6p plane, as in Fig. Q] The interior 
region of the ellipse violates m and hence represent 
gapped modes. Outside the ellipse the square root in 
is real, and is obviously satisfied as long as the left- 
hand side is negative. A negative left-hand side of la 
is represented by the region in-between the two branches 
of the hyperbola, p^ + Sp^ = M^ + SM"^ -I- A^, in Fig. ^ 
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For a positive left-hand side, which corresponds to the 
two regions on the left- and right-hand side of the hyper¬ 
bola, 03 can be squared, and four coupled inequalities 
linear in 5M and (5/r are obtained. These correspond to 
tangent lines of the ellipse. The hatched areas enclosed 
by the tangent lines, the hyperbola and the ellipse vio¬ 
late 113 and hence represent gapped modes. For each 
tangent line the intersection with the hyperbola coin¬ 
cides with the point on the ellipse. Inequality C3 is 
relevant if SM/p, ^ 1 — A^(/2 — M)“^/2, which is not 
the case for realistic values of the masses and chemical 
potentials. This is explicitely demonstrated by the ex¬ 
amples in Fig. 01 Inequality Ijl fill is therefore the relevant 
condition for gapless modes. In achieving this result, no 
further approximations than those leading up to were 
made. 

For the two-flavor color superconducting (2SC) phase, 
which is characterised by A^d ^ 0 and = Ads = 
= 0, one can use the fact that 6 M and M fi, 
so the gapless condition is approximately 


Aud < \SfJ.\ = -fiQ. (18) 


For the three-flavor color-flavor locked (CFL) phase, 
which is characterised by Ay ^ 0 and M ~ 5M ~ Ms/2, 
a series expansion of to first order in ^ yields 


Ais ^ |(5/i| + 


sign(J^) 

4 


Ml Ml 

8/2 p ’ 


(19) 


where 6 fj, = {fdia — /rs/3)/2. These well-known approxi¬ 
mative results are instructive at the qualitative level, but 
should not be used mechanically. See discussion below. 

Next we present a numerical example and therefore 
constrain the discussion to a specific parametrization 
of the model, as in 13 The momentum integral is 
regularized with a cut-off, A = 602.3 MeV. The cou¬ 
pling constants are = 2.319, Gd/Gs = rj = 0.75 

and K = 0. The current quark masses are m„ = 
nid = 5.5 MeV and rus = 112 MeV. Inserting these pa¬ 
rameters in the thermodynamic potential the gaps 
{(t)u,(t>dAs, Aud,Aus,Ads) can be determined by mini¬ 
mization of the free energy, while simultaneously neutral¬ 
izing all charge densities with fig, and /is. In Fig. [3 
the diquark gaps. Ay, and the gapless thresholds 
A® , are plotted vs. the temperature for a fixed value of 
the quark number chemical potential, /i = 500 MeV. In 
Fig.El the quark-quark quasiparticle dispersion relations 
are plotted for four different temperatures represented in 
Fig.El Observe that gapless dispersion relations exist iff 
Ay < A®.. Fig.H shows the graphical representation of 
the gapless condition for some quasiparticles represented 
in Fig.El see Fig.^for further information. 

We find that the difference between the approximative 
results and the exact gapless condition cni) is 

typically below 5% in the plane of temperature and quark 
number chemical potential. However, even a small error 


in A®^- could lead to qualitatively incorrect conclusions 
if Ay (T ~ 0) ~ A® , because Ay (T) are roughly con¬ 
stant at the low temperatures relevant for compact star 
evolution. The exact condition for gapless quasiparticle 
excitations presented here iHTni . which is the main result 
of this paper, is a safe alternative to the approximative 
results. Moreover, Fig. 0] is an accurate qualitative pic¬ 
ture of the prerequisites for gapless color-3 excitations in 
NJL models. 
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FIG. 1: Graphical representation of the gapless inequalities ITiH-ltTTt. For clarity only the first quadrant is shown. The third 
quadrant is a reflection of the first quadrant in the origin. In this figure an unreasonably large value of M has been used in 
order to emphasize the role of the tangent lines. Qualitatively, the gapped region can be represented by the rotated ellipse, see 
the text. The values of 5^i and 5M can be represented by a point in the 5M — Sfi plane. If this point is enclosed by the hatched 
area the dispersion relations are gapped. Otherwise a gapless dispersion relation exists. Here SM'^ = /i/(l + and 

= [A^ + (;t±M)2]l/2^ 



FIG. 2: Diquark gaps vs. the temperature at ^ = 500 MeV and rj = 0.75. AT is the threshold for gapless quasiparticle 
dispersion relations ITBll . i.e., gapless modes exist iff Aij < AT. The critical points where gapless db — sg, ub — sr and 
ug — dr quasiparticles appear are denoted with, respectively, A, B and C. The BCS result for the critical temperature of a 
superconducting condensate, T ~ 0.57A(r = 0), is indicated in the plot. This figure represents a cross-section of Fig. 5 in m. 
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FIG. 3: Dispersion relations at /i = 500 MeV and rj = 0.75 for four different temperatures. Gapless modes are denoted with 
bold lines. Gapped modes and modes of unpaired quarks are denoted with thin lines. Gompare with Fig.|^ 
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FIG. 4: Graphical representation of the gapless conditions (EJ-lO for some of the quasiparticle dispersion relations represented 
in Fig. 13 The values of 5M and Sfi are represented by bold points. If the center of a point is enclosed by the ellipse the 
corresponding quasiparticle have gapped dispersion relations, otherwise a gapless dispersion relation exists. 






















